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Preface 


This  report  logically  follows  previous  investigations  on  conditional  probabil- 
ity. The  first  report,  appearing  in  the  Monthly  Weather  Review  (July  1972),  was 
written  for  use  with  a single  predictor;  the  second,  appearing  as  AFSG  No.  354, 
provided  for  several  predictors  that  could  contribute  information  on  one  predictand 
by  linear  regression.  However,  for  two  predictands  occurring  jointly,  such  as 
the  ceiling  and  the  visibility  at  a target,  it  does  not  suffice  to  produce  individual 
predictions  for  the  conditional  probability  of  each.  To  meet  the  additional  require- 
ment, this  report  is  based  on  in-house  work,  done  at  the  request  of  Air  Weather 
Service  (RCS-7-11). 
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Conditional  Joint  Probabilities 


I.  INTRODUCTION 

The  problem  of  estimating  the  probability  of  concurrent  events  arises  fre- 
quently. One  of  the  most  persistent  problems  is  the  finding  of  probability  wherein 
both  ceiling  and  visibility  at  a station  jointly  will  have  values  above  their  thresh- 
olds. In  this  study,  the  same  predictors  will  serve  both  weather  elements,  although 
prior  ceilings  should  be  better  than  prior  visibility  for  the  prediction  of  subsequent 
ceiling,  and  prior  visibilities  should  serve  more  efficiently  for  subsequent  visibility. 

The  conditional  probability  of  two  jointly  occurring  predicted  events,  given  one 
or  more  predictors,  is  clearly  dependent  upon  the  conditional  probability  of  each 
individual  predictand.  It  is  also  dependent  on  the  correlation  between  predictands. 
Generally,  the  latter  correlation  is  conditioned  upon  selection  of  predictors. 


2.  FORMULA  FOR  CONDITIONAL  PROBABILITY 

Before  proceeding  further,  one  takes  for  granted  that  each  variable  has  been 
transformed  into  the  equivalent  normal  deviate  (END)  through  its  cumulative  fre- 
quency. The  variable  X is  transformed  into  the  END,  x(0,  1),  through  the  cumu- 
lative probability: 
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(1) 


P(^X)  P(>x)  — — f 


-s2  2 

e - 4 d? 


There  is  a convenient  approximate  formula1  for  determining  the  END  (x)  in 
terms  of  the  probability  of  exceedance,  [P(2:X)1,  as  follows: 

x k[t  - (a^  4 a^tj/d  + bjt  + b9t“)l  (2 


where 


a 2. 30753 
o 

at  0.27061 
b 0. 99229 
b„  0. 04481 


and  where 


k 1 , t 


I In  1/p  when  p - P(>X)  < 1/2 


k -1  , t = \fn  1/(1  - p)2 *  when  p = P(>X)  > 1/2 

If  yp  y2  are  the  predictand  END'S  and  Xj,  •••,  x the  predictor  END'S,  then 
the  problem  becomes  the  finding  of  conditional  probability: 

P(yl  - ylc-  y2  - y2c!xP  • • * - xn>  P(yl  " ylc-  y2  a y2c^ 

where  y^c  are  the  threshold  END-values  whose  joint  probability  of  exceedance 

is  to  be  found,  and 

X <xv  Xn>  . 

Correlation  coefficients  (CC's)  are  assumed  to  exist  between  each  pair  of 
variables,  with  each  CC  constant  throughout  the  distribution  of  the  variates. 


1.  National  Bureau  of  Standards  (1964)  Handbook  of  Mathematical  Functions 

Applied  Mathematics  Series,  55,  Government  Printing  Office,  Washington, 

D.  C.  20402,  pp  932-933. 


ii 


As  in  previous  work,  each  predictand  is  assumed  to  be  linearly  dependent 
on  'he  predictors.  Hence, 


yl  = allxl 


+ a,  x + b, p. 
In  n 11 


(3) 


y2  = an,x 


211 


a2nxn  + b2r)2 


where  rjp  17 2 are  normally  distributed  and  random  except  for  their  intercorrela- 
tion. The  a's  are  partial  regression  coefficients  and  the  b's  are  of  such  magnitude 
that  the  normality  of  Vp  y2  are  preserved. 

Introducing  the  vectors,  we  write 


‘11 


2 1 


^1 


X = 


In 


2n 


Equation  (3)  becomes 
T 


y!  ^1  * + hi'1! 


(4) 


y9  - .a,  x + b 9n 


2 '2 


Let  r . . be  the  CC  between  the  ith  predictor  and  the  jth  predictor.  Let  p...  be  the 


CC  between  the  ith  predictand  and  the  jth  predictor.  Then  the  coefficients  (A  p 
A2,  bp  b2>  can  be  obtained  in  terms  of  the  matrix: 


H, 


rll’ 

r12  r In 

r21’ 

r22  r2n 

rnl' 

\ 

rn2  * ’ rnn 

/ 

pll  \ 

— 2 = 

p In  i 

(5) 


(6) 


2.  Gringorten,  Irving  I.  (1976)  Multi-Predictor  Conditional  Probabilities.  AFSG 
No.  354,  AFSC  (AFGL),  FSAF,  Ilanscom  AFB,  24  pp. 
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In  terms  of  C,  Jip  _R.,,  we  have 


a2  = C'\ 


c'1  = c/|c| 


where  C is  the  adjoint  matrix  of  C,  and  |c|  is  the  determinant  of  C\  The  adjoint 
matrix  of  C is  a matrix  of  cofactors  (p.^)  as  follows: 


P.  . =(-Di+j  | My  | 

where  M--  is  the  submatrix  of  order  (n-1)  obtained  bv  deleting  the  ifb  row  and  the 
i] 

jth  column  of  jC.  Also 


i -a/  rx 


1 - a2t  r2 


Hence,  from  Eq.  (4) 


nl  = <yl  ‘ allxl  ' alnxn>  bl 
12  = (y2  - a21xl  •••  ' a2nxn>  b2 


Clearly, 


P(yl  £ ylc-  y2  -y2cl^)  p(,,l  - ^lc’  n2  - 11 2c'  • 


Since  r)2  have  Gaussian  distributions,  their  joint  probability  is  dependent  upon 
the  CC  ( p ) between  them  in 


9 2 

£ “ - 2p£q  + o 


p(hl  s nlc.  n2  > n2c) 


2*  \ l - p nlc  ')2c 


J ^ e p ' d£  dp  (12) 


To  find  p,  one  needs  to  recognize  it  as  the  expected  value  of  the  product  of 
r/j  and  r\2-  From  Eq.  (11): 


blb2,1l'12  >V’2  - ~ alixrv2  - E a2ixhVl  * E L alixix,a2.i  ’ 


bjbgp  po  - AjT  R2  - A,1  iit  +A^  CA2  (14) 

where  o i.s  the  unconditional  CC  between  v,  and  v0.  Finally, 

p --  (po  - AjT  CA2)/b1b2  (15) 

or,  in  purely  algebraic  terms 

p - £ a , • a „ r . . 
po  . 1 i 2 1 1 1 

p - M-  - ■ <ii 

V1  alipli  V1  ■ l a2ip2i 

Equation  (15)  or  (16)  can  be  recognized  as  the  partial  correlation  between  y^  and 

y„,  a result  that  is  a consequence  of  the  absence  of  observational  x -values  in  the 
4 3 

computation  of  p . In  this  report,  therefore,  the  terms  conditional  C'C  and 
partial  CC  are  interchangeable. 

3.  SAMPLE  STUDIES 


3.1  Ceiling  hikI  Visibility  in  Southern  Germany 

This  case  concerns  probability  at  Nurnberg  (49°30'N,  11°5'E),  given  observa- 
tions at  Kitzingen  (49°42'N,  10°6'E).  The  distance  between  stations  is  74.21  km. 

4 

In  Model  B,  the  CC's  between  ceiling  or  visibility  at  one  station  and  ceiling  or 
visibility  at  another  station  are  given  by  the  formula: 


, • — [sir 

O 7T  L 


<1  - a"  - a V 1 


3.  hawrance,  A.I.  (1976)  On  conditional  and  partial  correlation.  The  American 

Statistician  30:146-149. 

4.  Oringorten,  Irving  I . ( 1978)  Modelling  for  the  Probability  Distribution  of  Ceil 

ing  and  Visibility  at  a Point,  Along  a Line  of  Travel,  or  in  an  Area,  in  the 
Southern  Region  of  West  Germany,  Unpublished  Report,  submitted  to  Hq. 
AWS,  26  April  1973. 


where  s'  is  the  distance  between  stations,  r is  the  parameter  known  as  the  scale 
distance,  and  p is  a correction  factor,  mainly  due  to  observational  errors,  equal 
to  0.  85  for  ceiling  versus  ceiling,  and  0.  85  for  visibility  versus  visibility.  For 
ceiling  versus  visibility,  the  value  of  po  is  0.  17. 

For  the  German  mountain  stations,  r - 6.  24  km  in  winter.  Hence,  for 
s'  74.  21  km,  Eq.  (17)  gives  the  January  CC's  that  compose  the  matrix  C and 
vectors  RR  and  Ry  for  ceiling  height  (H)  and  visibility  (V)  as  follows: 


V 75  ' 

l o-i 5 1 

From  Eqs.  (7)  to  (11),  we  find 
| C 1 = 0.9711 


0.  85  0.  17 


, 0.  1 7 0.  85  , 


— H 


/ 0.  85 

-0.  17  ^ 

C*1 

( 1.  2255 

-0.  2451  \ 

1-0.  17 

0.  85  j 

y -0.  2451 

1.2255  ) 

( 0.  8824 

A,, 

° ) 

( ° J 

— V 

0.  8824  1 

bjj  - 0.  582  , by  = 0.  582 


Finally, 


H = (y  - 0.  8824  xR)/0.  582 

(18a) 

v = (y  - 0.  8824  Xy)/0.  582 

(18b) 

Fortuitously,  ceiling  depends  on  ceiling  predictor  only;  visibility  depends  on  visi- 
bility predictor  only.  From  Eq.  (15)  or  Eq.  (16) 


p - 0.  118 


To  find  the  joint  probability,  Eqs.  (18a  and  b)  must  be  used;  the  values  of  xR 
and  Xy  corresponding  to  the  observations  of  ceiling  and  visibility  at  Kitzingen 
must  be  known,  and  the  threshold  values  of  yR  and  y at  Nurnberg  be  selected. 
For  this  report,  the  END'S  are  found  through  the  cumulative  frequencies  of  the 


12 


r 


-v—iy 


ceilings  and  the  visibilities  at  the  two  stations,  Kitzingen  and  Nurnberg  (Figures  1 
and  2). 

Suppose  that  the  ceiling  and  visibility  at  the  predictor  station,  Kitzingen,  are 
known  to  be  300  meters  and  1 km,  respectively.  Then,  from  Figures  1 and  2, 

Xjj  -1.00  and  -1.60.  Therefore  Eq.  (18)  becomes 

1H  1.719  yH  + 1.339 

r)v  = 1.  719  yy  + 2.  427 

Table  1 gives  values  of  y„  and  y,r  for  several  threshold  values  of  ceiling  and  visi- 
bility at  Nurnberg,  the  predictand  station.  The  entries  for  joint  conditional  prob- 
ability were  found  by  numerical  solution  of  Eq.  (12)  on  a Hewlett-Packard  9810A. 

3.2  Joint  Conditional  Probability  of  24-hour  Rainfall  in  \i*h  hn^land 

The  selected  example  is  that  of  the  joint  probability  of  24-hour  rainfall  at 
Boston,  Massachusetts,  and  Portland,  Maine.  The  question  concerns  whether  or 
not  there  has  been  rain,  and  how  much,  at  the  more  westerly  stations:  Providence, 
Hartford,  and  Burlington.  The  five  stations  are  located  as  follows: 


Boston 

Portland 

Providence 

Hartford 

Burlington 


42°22'N  71°02'\V 

43°39'N  70°19'U 

4 1°44'N  71°2fi'\V 

41°56'N  72°41'\V 

44°28'N  73°09'\V 


From  the  formula  for  distance  s'  (kilometers)  between  the  it h and  it h stations: 


111.12  b 


(b  in  degrees  and  tenths) 


where 


cos  b - sin  d-  sin  d + cos  o-  cos  d cos  (X.  - X.) 

1 .1  1 .1  i .1 

and  where  d represents  latitude  and  X longitude.  For  distances  between  stations, 
see  Table  2. 

3.2.1  JANFARY 

Previous  records  gave  the  probability  distribution  of  24-hour  precipitation  in 
New  England  (January),  producing  the  curve  of  rainfall  versus  END  (f  igure  3). 
The  26-year  record  (1952-1977)  yielded  an  average  frequency  of  0,3801  for 
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KITZINGEN 


NURNBERG 


NURNBERG 


KITZINGEN 


0001 


CEILING  HEIGHT  (METERS) 


Figure  1.  Cumulative  Frequency  of  Ceiling  Height 
(January),  0900  I at  Kitzingen  and  N'urnberg 


NURNBERG 


* TZlNGEN 


0 2 5 8 1012  15  20  25  30  35 

VISIBILIT r ikM) 


Figure  2.  Cumulative  Frequency  of  Visibility 
(January),  0900  L,  at  Kitzingen  and  Xurnherg 


Table  1.  Threshold  Ceiling  and  Visibility.  Correspondmg  values  (yH,  yv),  un 
conditional  probabilities  [P(>yH),  P(>yy)l  , resulting  values  (pjj,  n\r)  for  the 
specific  predictors  at  Kitzingen  (ceiling  300  m,  visibility  1 km),  and  the  condi 
tional  probabilities  [P(ar)jj),  P(>0y)]  (Estimates  are  shown  of  the  joint  condi- 
tional probabilities  when  the  conditional  correlation  coefficient  is  0.  118) 


Visibility  >1  km  3:5  km 

yy  -1.  50  0.31 

P(>yv)  0.  933  0.378 

rj„  -0.15  2.96 


Ceiling 

lit 

P(SvH) 

nH 

P(^H> 

p(>nv) 

0.  560 

> 1 50  m 

-1.  40 

0.  92 

-1. 07 

0.  86 

0.  489 

>300  m 

-0.  72 

0.  76 

0.  10 

0.  46 

0.  275 

>600  m 

-0.  18 

0.  57 

1.03 

0.  15 

0.  095 

> 1 500  m 

0.35 

0.  36 

1. 94 

0.  026 

0.  017 

>9000 m 

0.  94 

0.  17 

2.  95 

0.  0016 

. . . 

0.  0012 

0.  0010 


0. 000084 


Table  2.  Distances  Between  New  England  Stations 
(in  kilometers) 


Station 

Boston 

i 

Portland 

Providence 

Portland 

Providence 

Hartford 

Burlington 


77.  8 
144.  3 
289.  1 
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Figure  3.  Transformation  of  New  England  24-hour  Precipitation  (January) 


rainfall  >0.25  mm.  The  frequencies  and  joint  frequencies  of  rain  >0.25  mm  were 
as  shown  (Table  3).  The  distances  between  stations,  coupled  with  a previously 
determined0  scale  distance  of  r 9.45  km  for  24-hour  precipitation  gave  COs  be- 
tween stations  (Table  4).  These,  in  turn,  together  with  the  marginal  frequencies 
(Table  3),  as  used  in  Eq.  (12),  gave  the  estimates  of  joint  probabilities.  On  the 
whole,  the  estimates  are  low  by  an  average  of  0.02,  suggesting  that  the  parameter 
size  (r  = 9.45  km)  could  have  been  chosen  somewhat  larger,  at  least  for  winter 
precipitation  equal  to  or  greater  than  0.  25  mm. 

For  precipitation  >21.  84  mm,  the  marginal  and  joint  frequencies  were  as 
shown  (Table  5).  If,  again,  the  CC's,  as  estimated  from  Model  B with  r 9.45  km, 
are  used  in  Eq.  (12)  the  estimates  of  joint  probability  become  the  figures  in  paren- 
theses (Table  5).  This  time  the  average  bias  is  less  than  0.001. 

3.2.2  JULY 

Previous  records  gave  the  probability  distribution  of  24-hour  rainfall  in  July 
as  shown  (Figure  4).  The  average  frequency  is  0.305  of  rainfall  >0.25  mm.  The 
marginal  and  joint  frequencies  >0.25  mm  were  as  shown  (Table  6). 

5.  Gringorten,  Irving  I.  (197(1)  Areal  Coverage  Estimates  by  Stochastic  Modelling, 
ERPNo.  573,  AFSC  (AFGE),  t'SAF,  Hans  com  AFB,  50  pp. 


If, 


0 


i 


Station 

Boston 

Portland 

Providence 

Hartford 

Boston 

0.  385 

0.300 

0.  328 

0.  318 

0.  284 

(0.28) 

(0.31) 

(0.  29) 

(0.  28) 

Portland 

. . . 

0.  3 54 

0.  279 

0.  275 

0.  284 

(0.25) 

(0.  24) 

(0.  27) 

Providence 

... 

0.  367 

0.  315 

0.  259 

(0.  29) 

(0.  26) 

Hartford 

0.364 

0.  273 
(0.  27) 

Burlington 

. . . 

. . . 

0.  438 

Numbers  in  parentheses  are  Model  B probability  estimates. 


Table  4.  Model  B Estimates  of  Correlation  Coefficients  of  24-hour  January 
Precipitation  Between  Pairs  of  Stations 


Table  5.  Marginal  and  Joint  Frequencies  of  24-hour  January  Precipitation 
221.  84  mm.  (Single-station  frequencies  are  shown  on  the  diagonal) 


Station 

Boston 

Portland 

Providence 

Hartford 

Burlington 

Boston 

0.0347 

0.  0099 
(0.  Oil) 

0.  0236 
(0.020) 

0. 0136 
(0.  Oil) 

0. 00124 
(0.  0020) 

Portland 

0.  0298 

0.  0136 
(0.  012) 

0. 00868 
(0.  0068) 

0. 00124 
(0.  00248) 

Providence 

. . . 

0.  0422 

0. 0136 
(0.  014) 

0. 00248 
(0.0019) 

Hartford 

. . . 

... 

. . . 

0.0223 

0. 00248 
(0.  0016) 

Burlington 

... 

0.0037 

Numbers  in  parentheses  are  Model  B estimates. 


Figure  4.  Transformation  of  New  England  24-hour  Rainfall  (July) 


The  previously  determined''  scale  distance  for  July  24 -hour  rainfall  is 
r = 5.  19  km.  Coupled  with  the  distances  between  stations  as  determined  above, 
the  CC's  between  stations  are  as  estimated  (Table  7).  The  CC's,  in  turn,  together 
with  the  marginal  frequencies  in  Table  (i,  when  used  in  Eq.  (12),  yielded  the  joint 
probabilities  in  parentheses  (Table  !>).  The  latter  show  a bias  of  only  0.003  or 
3/10  of  1 percent. 

For  rainfall  ?27.  94  mm  (END  2-2.0),  the  marginal  and  joint  frequencies  were 
as  shown  (Table  8).  The  Model  13  estimates  of  joint  probability  are  the  bracketed 
figures.  They  show  a bias  of  -0.001,  and  an  RMSE  of  0.001. 

This  supports  our  reliance  upon  Model  13.  The  scale  distances  of  Model  B 
were  used  to  estimate  the  CC's  that  were  vised  in  Eqs.  (4)  to  (15)  in  order  to  find 
the  partial  regression  coefficients  and  other  parameters  (bj,  - b9),  and  finally  the 
conditional  probabilities  (Section  3.2.3). 
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Table  7.  Model  B Estimates  of  Correlation  Coefficients  of  24-hour  July  Rainfall 


St  at  ion 

Boston 

Portland 

Providence 

Hartford 

Burlington 

Boston 

1. 000 

0.  707 

0.  851 

0.  725 

0.  4 04 

Portland 

1. 000 

0.  5f>5 

0.  494 

0.  543 

Providence 

1. 000 

0.  798 

0.387 

Hartford 

1. 000 

0.  472 

Burlington 

... 

1. 000 

Table  8.  Marginal  and  Joint  Frequencies  of  24-hour  July  Rainfall  2-27.  94  mm 


Station 

Boston 

Portland 

Providence 

Hartfo  rd 

Burlington 

Boston 

0.  0174 

0.  000775 
(0.  0031) 

0. 00574 
(0.  0004) 

0.  00326 
(0.  0047) 

(0.  0015) 

Portland 

0.  0149 

0.  00078 
(0.0027) 

0. 00326 
(0.0023) 

0.  000775 
(0.0017) 

Providence 

... 

0.0223 

0.00574 
(0.  0009) 

(0.  0014) 

Hartford 

. . . 

CT 

c'i 

o 

© 

0.  00202 
(0.  0019) 

Burlington 

— 

0. 0136 

Numbers  in  parentheses  are  Model  B estimates. 
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Table  9.  Marginal  and  Joint  Conditional  Probabilities  of  24-hour  Precipitation 
at  Boston  and  Portland  (January) 


Actual  precipitation:  None  at  Boston,  0.25  mm  at  Portland 


On  3 Jan 

1953,  x1 

= 2. 24  x2  1. 

97.  *3 

1.  51* 

>0.  25 

0.4 

-4.  46 

0. 9999959 

-2.  79 

0.  9974 

0. 

9974 

>4.  2 

1.0 

-2.  85 

0.  9978 

-1.  72 

0.  957 

0. 

955 

£21.  8 

2.  0 

-0.  17 

0.  57 

0.  07 

0.47 

0. 

342 

Actual  precipitation:  30.  5 mm  at  Boston  and  at  Portland 


See  text 


Table  10.  Marginal  and  Joint  Conditional  Probabil it ies  of  24-hour  Rainfall 
at  Boston  and  Portland  (July) 


For  details  on  the  values  of  x^,  x^,  x^,  see  discussion  in  text. 
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3.  2.  3 JANUARY  PARAMETERS  AND  CONDITIONAL 
PROBABILITIES 

Using  B,  P as  subscripts  for  Boston  and  Portland,  respectively,  and  1,  2,  3 for 
Providence,  Hartford,  and  Burlington  respectively,  one  finds  that  the  solutions  of 
Eqs.  (5)  to  (16)  produce  the  following  results: 


1 r12  r13 

r2 1 1 r23 

r3 1 r32  1 


0.889  0.652 


0.  889  1 0.  704 

0.652  0.704  1 


0.  918 


0.  717 


from  which 


C =0.  1051 


0.5044  -0.4300  -0.0261 

-0.4300  0.5749  -0.1244 

-0.0261  -0.1244  0.2097 


4.799  -4.091  -0.248 

-4.091  5.470  -1.184 


-0.248  -1.184  1.995 


Ab  I 0.  0566 


0. 0611 


bB  = 0.  373 


bp  0.  561 


From  Eqs.  (11)  the  END'S  become 


hB  = <yB  * °-  763x1  - 0.  0566x2  - 0.  161x3>/0.  373 
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p = <yp  - 0.  520x1  - 0.  0611x2  - 0.  449x3)/0.  561 


From  Eq.  (16) 

p = (0.  838  - 0.739)/(0.373)(0.  561)  = 0.473 

which  makes  the  conditional  CC  substantially  reduced  from  the  unconditional  of 
0.  838. 

When  q r )_  are  found,  from  the  above  equations,  for  several  END'S  of 

threshold  rainfall  (yB,  yp),  for  specific  END’S  of  predictor  information  (Xj,  x2, 
x3>,  the  results  are  as  shown  (Table  9).  On  1 January  1952,  there  was  no  rain 
at  Providence,  Hartford  or  Burlington.  Since  the  probability  of  no  rain  at  a single 
station  is  0.  66,  the  END  was  chosen  for  half  this  probability  (that  is,  Xj  = x2  = 
x,j  = -0.44).  On  3 January  1953,  there  was  27.9  mm  of  precipitation  at  Providence, 
24.  1 mm  at  Hartford,  11.  9 mm  at  Burlington  (which  Figure  3 converts  to  x:  = 2.  24, 
x2  = 1.  97,  x3  = 1.  51). 

3.  2.4  JULY  PARAMETERS  AND  CONDITIONAL  PROBABILITIES 

In  similar  procedure,  the  calculations  give: 

nB  = (yR  - 0.  746  xt  - 0.  060  x2  - 0.  147  Xg)/0.  503 

np  = (y  - 0.  459  Xl  + 0.  058  x,  - 0.  393  x3)/0.  746 


and 


p = 0.  469 

Again  the  conditional  CC  is  reduced  but  still  significant. 

On  1 July  1952,  there  was  no  rain  at  the  three  predictor  stations.  Since  the 
probability  of  no  rain  is  0.  69,  the  END  is  taken  at  one  half  the  probability.  Thus 

X1  = x2  = x3  " -°-40- 

On  2 July  1953  there  were  0.  76  mm  of  rain  at  Providence,  1. 5 mm  at  Hartford, 
and  5.  1 mm  at  Burlington.  Hence  Xj  = 0.  70,  x2  = 0.  81,  x3  = 1.  09. 

Table  10  shows  the  results  on  the  conditional  probabilities,  similar  to  those  in 
Table  9. 


4.  CONCLUSIONS 


The  conditional  joint  probability  of  two  events  that  are  conditional  upon  one  or 
more  pieces  of  predictor  information  can  be  related  to  the  conditional  probabilities 
of  each  event  separately,  given  the  conditional  CC  between  them.  To  find  the 
individual  conditional  probability  of  each  predictand,  the  corresponding  equivalent 
normal  deviate  is  determined,  using  an  equation  involving  the  deviates  of  predic- 
tand and  predictors.  The  conditional  CC  in  this  work  is  synonymous  with  partial 
correlation  coefficient;  it  is  found  in  terms  of  the  unconditional  CC  corrected  by 
terms  involving  CC's  between  predictand  and  predictors,  and  also  between  pre- 
dictors and  predictors. 

While  the  conditional  CC's  have  proved  smaller  in  the  examples  cited,  it  is 
conceivable  that  they  could  have  remained  unchanged.  Consider  the  first  example, 
ceiling  and  visibility  in  the  mountain  regions  of  Southern  Germany.  The  effect  on 
the  joint  conditional  probability  has  been  that  of  estimating  only  a slight  improve- 
ment in  the  joint  probability  over  the  product  of  the  two  separate  probabilities. 

These  estimates  cannot  be  verified  by  sampling. 

In  the  case  of  New  England  rainfall,  the  viability  of  Model  B was  first  estab- 
lished by  its  estimates  of  joint  unconditional  probabilities.  The  Model  B estimates 
of  the  CC's  produced  substantially  reduced  conditional  CC's  between  the  predictands. 
The  resulting  estimates  of  joint  conditional  probability,  however,  in  some  instances 
were  substantially  higher  than  the  product  of  the  separate  conditional  probabilities. 
The  most  significant  result  can  be  seen  in  the  wide  variations  in  estimates  of  con- 
ditional probabilities  both  separate  and  joint,  depending  on  the  predictor  values; 
that  is,  depending  on  what  the  observers  see  at  the  predictor  stations. 
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